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ABSTRACT 

We present an analysis of HCO+ (J = 3 -»■ 2) and H 13 CO+ (J = 3 — >• 2) 
observations of the massive (M~ 20 M ) submm/IR source IRAS 01202+6133 
located on the periphery of the H II region KR 120 (Sh 2-187). The HCO+ line 
profile has a classic blue-asymmetric shape with the optically thin H 13 CO + line 
peaking at the position expected if the HCO + line arises from a combination of 
self-absorption and infall motion. We have modified existing analytic radiative 
transfer models to allow for the fitting of submm/mm line profiles that have both 
self-absorption features and optically thin wings and applied these models to our 
HCO+ spectrum of IRAS 01202+6133. We conclude that it is a young Class I 
YSO with a substantial envelope undergoing slow infall and having some outflow 
motions. The young age of the H II region rules out a "collect and collapse" sce- 
nario. While we cannot eliminate the possibility that IRAS 01202+6133 formed 
spontaneously at its current location, considering its early evolutionary state and 
its proximity to the H II region we think that the formation of IRAS 01202+6133 
was triggered by the expansion of KR 120 (Sh 2-187). 

Subject headings: ISM: molecules - ISM: kinematics and dynamics H II regions 
- stars: formation - Line: profiles 



Introduction 



The molecular gas surrounding an H II region is thought to be a place where star 
formation can be induced. Such triggered s tar formation can aris e from the overpres- 



surization of existing density enhance ments (Thompson et al. 2004) or throug h the col- 



lapse of a swept up layers of material (jElmegreen 



1998; 



Deharveng et al.l 120031 ) . This pa- 



per presents an analysis of submm spectroscopic observations of the Class I Young Stel- 
lar Object (YSO) IRAS 01202+6133 which is located on the periphery of the relatively 
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nearby (d = 1.44 ± 0.26 kpc) 


H II 


Kalas & Reich 


1980; 


Sharpless 


1959) 



sive (M = 21 ± 9 M ), luminous (L ~ 5600 L ) YSO to the H II region makes it a strong 
candidate for an example of triggered star formation. The submm spectra presented here 
strengthen this scenario as they clearly show infall is occurring as would be expected for a 
very young YSO. We believe that IRAS 01202+6133 is one of the clearest examples to date 
of an H II region triggering the formation of another massive star. 



Our o bservations are des c ribed in §HJ In §[3]we review the analytic models of lMyers et al. 



( 119961 ) and iDeVries fc Myersl ( 120051 ) and develop new varieties of these models that allow for 
both infall and outflow motions. Physical properties of the YSO derived using these models 
are presented and discussed in § 0] followed by our conclusions in § 



2. Observations 

In order to explore t he kinematics of I RAS 01202+6133 we obtained observations of 



the commonly used (e.g., ISun &: Gaol 120091 ) optically thick infall tracer molecule HCO + 
(J = 3 -> 2; 267.6 GHz) and its optically thin isotopologue H 13 CO+ (J = 3 ->• 2; 260.3 
GHz) at the James Clerk Maxwell Telescope (JCMT) in the 06B semester. Single pointing 
observations (beamsize ~ 20" FWHM) were made using the A3 receiver and ACSIS backend 
with a frequency resolution of 61 kHz (0.068 km s -1 ) and an integration time of ~ 40 minutes 
for each line. The reduced spectra are shown in Figure [2j Both have a noise level of ~ 0.1 K 
resulting in a peak S/N~ 30 for the HCO + spectrum and ~ 3 for the weaker H 13 CO + 
spectrum. 

If there was no infall motion associated with IRAS 01202+6133 we would expect the 
HCO + spectrum to be a single self-absorbed profile with sym metric red and blu e peaks 



around a local minimum at some systematic velocity. Following iMyers et al.l (119961 ) . in the 
case of where there is infall there will be density and excitation temperature (T ex ) differences 
between the central regions and the outer parts of the molecular core. Blueshifted emission 
from the rear of the core is not strongly absorbed because it passes through the central (high 
T ex ) region then through a highly Doppler shifted front portion of the core. In contrast the 
redshifted emission from the front of the core is absorbed by surrounding gas that is at both a 
comparable velocity and similar (low) T ex . The overall effect of the infall motion is to create 
a "blue-asymmetric" line profile where the self-absorption line is asymmetric with a stronger 
blue peak. Such a line profile could also form from multiple unrelated emission components 
being aligned along the line of sight or from rotational motion, but in both cases we would 
expect to see a similar profile in the optically thin line. The fact that our H 13 CO + spectrum 
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peaks up at the same velocity as the minimum in the HCO + spectrum makes us confident 
that we are looking at a single self-absorbed line profile with the asymmetry arising from 
infall motions. 



3. Analytic Radiative Transfer Modeling 



Myers et al.l ( 119961 ) and lDeVries &: Myersl (120051 ) developed basic one- dimensional (ID) 



analytic radiative transfer models to derive kinematic data from mm/submm spectra of 
star-forming molecular cores. In this section we first review these models and then present 
modified versions that better model cases where there is simultaneous infall and outflow 
motions. 

All of the models are based on the general solution to the ID radiative transfer equation 
for a homogeneous medium with total optical depth r , 

/•TO 

T B = T ie - T0 + / J(T)e- T dr (1) 
Jo 

where T B is the brightness temperature of the outgoing radiation, Tj is the brightness tem- 
perature of any incoming radiation, J(T) = T /[exp(T /T) — 1] and T = hu/k, where v is 
the frequency of the transition and T is the temperature. 

If J(T) is a linear function of r then equation [1] can be integrated to obtain, 

T B = T ie - T0 + (J 2 - Ji) 1 " 6 T ° + Jx - J 2 e- T0 (2) 

TO 

where J% and J 2 are the values of J(T) at the two ends of the path with optical depth To- 
All of the models presented in this section are based on equation [2] with variations due to 
how J(T) varies with r in detail. 



3.1. Two- layer model 



The two-layer model of Myers et al.l (119961 1 has "front" and "rear" c onstant T PX layers 
movi ng toward each other each at an infall velocity, Vi n (see Figure 1 of iDeVries fc Myers 



20051 ). 



The line brightness temperature of the two-layer model is 

AT B = J(7»[l - e' T f\ + J(T r )[l - e - Tr \e- T t - J{T b )[l - e~ T ^ 



(3) 
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where Tf and T r are the front and rear excitation temperatures and T b is the background 
temperature. The optical depths of each layer, 77 and r r are given by 



t exp [-(v - v LSR - v in ) 2 /2a 2 ], 
t exp [-(v - v LS r + v in ) 2 /2a 2 } 



(4a) 
(4b) 



where r is the line center optical depth, v lsr is a average line-of-sight velocity of the system, 
and a is the velocity dispersion of the observed molecule. 

For a given choice of T b (usually T b = 2.7 K) the two-layer model has six adjustable 
parameters: Tf, T r , r , Vlsr, and a. By making some assumptions about the expected 



density and temperature structure in an infalling core iMyers et al.l (119961 ) was able to reduce 
the number of parameters in the two-layer model to five, replacing T r and Tf by a kinetic 
temperature, T^. This version of the model also used the non-thermal velocity dispersion 
as an input, g^t = (c 2 — kTk/rriobs) 1 ^ 2 where m b s is the mass of the observed molecule. 
H ere we presen t the formulas used to calculate the model parameters and refer the reader 



to 



Myers et al.l (119961 ) for their derivation. 



The excitation temperature of each layer is determined using the following equations, 



Tj_ T b + ( y 4:T /f3)(n f )/n max 
T k T k + (4T / 75) (n/) /n max ' 
Tr_ T b + (4T /(3)(n r )/n max 
T k T k + {4T /P)(n r )/n 

max 

where (5 = [1 - exp(-r )]/r , and (n f )/n max and (n r )/n max are given by, 



(5a) 
(5b) 



(n f ) 



n, 



(n r ) 



'1 - 



6 



6 6 



T 



n, 



r o r o T o 



(6a) 
(6b) 



3.2. Hill model 



In the "hill" model of iDeVries fc Myerd (120051 ) J(T) increases linearly from the Tq to 



T p over a front optical depth (7/), and then decreases from T p to T over a rear optical 
depth {T r ). The shape of a J{T) plot against r gives the model its name (see Figure 1 of 
DeVries fc Myersll2005l ). As with the two-layer model the front and rear sections each have a 
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line center optical depth (to), are moving toward each other with an infall velocity V{ n , and 
have a velocity dispersion a. 

The line brightness temperature of the hill mode0 is given by, 



AT B = [J{T P ) - J (T )] 



1 - e' T f e~ T f ( 1 - e~ Tr 



Tr- 



+ [J(T ) - J(T b )] (1 - e- T ^) . (7) 



For a given choice of T&, the hill model has six adjustable parameters: T , T p , r , vlsr, 
and a. Equations Hal and libl are used to determine r/ and r r . 



3.3. Analytical models with outflow 

While the two-layer and the hill models are able to reproduce the central portion of the 
observed blue-asymmetric line profile they are not able to model the extended wings of the 
line profile observed beyond ~ ±2 km s _1 (see Figure [3]). In orde r to correctly model the 
entire line we follow the approach used in the iMyers et al.l (119961 ) analysis of L 1527 and 
incorporate a central outflow region into both the two-layer and hill mod els. The resulting 
varia tion of T ex with r is shown schematically in Figure H] (cf. Figure 1 of iDeVries fc Myers 
2005h . 



3.3.1. Two-layer model with central outflow (2L-0 model) 

With the addition of a central outflow region to the two-layer model the line brightness 
temperature (equation |3]) becomes, 

AT B = J(T f )[l-e- T f] + J(T out )[l-e- T ^]e- T f + J(T r )[l - e^e^f^ (8) 
-J(T b )[l -e~ T f~ Tou 



out 1 r 



where T out and T out are the optical depth and excitation temperature of the outflow. 

The only additional parameters needed in the 2L-0 model are To y0 ut and a ou t the line 
center optical depth and the velocity dispersion of the outflow region. Both are used to 
calculate T out using, 

Tout = T ,out exp [-(« - v LSR ) 2 /2al ut ] . (9) 



1 The equation in DeVries fc Mver si ( 2005) is typeset incorr ectly in the Astrophysical Journal. The correct 
version is found in the arXiv version |DeVries k Mversll2004l ) 
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In our model we set T out equal to the average of T r and Tj calculated using T out fl in equa- 
tions [5a] [5bl EH and [6b] as appropriate. 

As a test of our version of the 2L-0 model we were abl e to exactly match the model 
H 2 CO (J = 2 — > 1) profile of L 1527 shown in Figure 3 of iMyers et al.l (119961 ) using the 
outflow and model parameters given in the figure caption. 



3.3.2. Hill model with central outflow (Hill-0 model) 
The brightness temperature of the hill model with a central outflow is 



AT B = (J p - J ) 



1 - e- T f ( 1 - e -Tr) e -Tf-Tout 



+ (J - J 6 )(l-e- T '- To » t - 7V ) (10) 



where r out is the optical depth of outflow. Since we set the excitation temperature of the 
outflow equal to T p (see Figure HJ) the only additional parameters are r 0tOUt and o out . As in 
the 2L-0 model, r out is calculated using equation EB 



4. Results & Discussion 

4.1. YSO Properties 

We applied the 2L-0 and the Hill-0 models to our observed HCO + spectrum (see Figures 
[5] and [6] respectively) . The best-fit parameters are given in Tables [1] and |2j Uncertainties 
were derived using a Monte Carlo method; an initial best fit of the model was obtained, 
then random noise at the same level as the observed line profile was added to the best fit 
and the resulting simulated line profile was then fit generating a new set of parameters. The 
standard deviation of spectral line parameters after 1000 repeats was combined with the 
mean fit uncertainties to produce the final error estimates. 



As noted by IMyers et al.l (119961 ) the two-layer model tends to output profiles that have 
narrower peaks and flatter central troughs than observations, and both of these characteris- 
tics can be seen in the 2L-0 best-fit model shown in Figure [5j In contrast the Hill-0 model 
does a much better job in fitting the observed peak widths and the structure of the central 
absorption trough. 



DeVries &: Myersl (120051 ) showed that Hill style models were able to obtain accurate 
estimates of infall velocities (rms error 0.01 km s -1 ) from two-peak profiles in contrast 
with two-layer models, which tended to underestimate the infall velocity by factors of ~ 2. 
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Consistent with this finding, we see that the best-fit 2L-0 Vi n is a factor of 3.5 lower than the 
Hill-0 value. While we adopt the numerical results of the Hill-0 model for IRAS 01202+6133 
we note that, regardless of the exact numeric values, both the 2L-0 and Hill-0 results 
are consistent with an object with a substantial envelope (r ~ 10) undergoing slow infall 
{V in < or)- 



4.2. Formation Scenarios 



There are two basic scenarios for the formation of IRAS 01202+6133, either its formation 
was triggered by the action of the expanding H II region or it formed spontaneously and its 
proximity to the H II region is coincidental. 

To gauge the likelihood of the first scenario we can first look at the various time scales 
associated with the YSO-H II regio n system. The canonic al duration of the Class I stage of 
YSO evolution is of order 10 5 years JWard-Thompsonll2002h . The fact that IRAS 01202+6133 
still has a substantial envelope and exhibits slow infall suggests that it is still a fairly young 
Class I YSO. The presence of outflow motions does not constrain the ag e as they are know n 
to be associated with even the very earliest stages of YSO evolution f lAndre et al.lll993l ). 
Assuming that any Class stage would be very rapid (~ 10 4 years), w ^e conclude that the 
age of the YSO is Tyso < 10 5 years. KR 120 is a fairly young H II region, IJoncas et al.l (119921 ) 
estimate an age t hii « 1 - 2 x 10 5 years based on kinematic and photodissociation region 
models for the region. So it appears that a necessary condition for triggered formation, that 
t yso < thii, is just satisfied in this case. 

As mentioned in § (TJ triggered star formation around H II regions can arise from the 
overpressurization of existing cores or through the collapse of a swept up layers of mate- 
rial. To investigate th e latter "collect and collapse" scenario we use the equation from 
Whitworth et al.l ( 119941 ) giving the timescale at which fragmenation and collapse will occur 
in a swept-up shell around an expanding H II region: 



t 



fragment 



1.56 Myr 



a s 


7/11 


r Qo 


-l/n 


n 


-5/11 


0.2 km s" 1 . 




_10 49 ph s~\ 




LlOOO cm" 3 . 





where a s is sound speed of the molecular cloud, Qo is the H I ionizin g flux, and no is the 
densi ty of the molecular cloud. The exciting star of KR 120 i s BO V (lArvidsson fc Kerton 
201 ll ) corresponding to Qo=10 47 ' 4 p h s -1 using the ICrowtherl (120051 ) calibration. For n we 



fll992h 



use the value n = 1000 cm -3 from IJoncas et al 
is approximately between 2 Myr with sound speed a 
clouds) and up to a 6 Myr upper limit with with a s - 



A fragmetation collapse timescale 
= 0.2 km s -1 (very cold molecular 
1 km s _1 (warm molecular /neutral 
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clouds). Since this fragmentation timescale is 10 times longer than thii we can eliminate 
the "collect and collapse" triggered star formation scenario. 

The morphology of the molecular cloud surrounding KR 120 (see Figure [7]) is consistent 
with the alternative scenario of triggered collapse of a pre-exisiting molecular core. The 
highest column density gas (e.g., Tb > 50 K km s -1 ), where one would expect to find 
molecular cloud cores, lies in a ridge to one side of the H II region. In this picture a shock 
front being driven into this gas by the expanding H II reg ion has overrun and compressed 



an p re-existing core initiating the star formation process (IBertoldil Il989l ; iThompson et al. 



20041 ). 



We cannot rule out a spontaneous star formation scenario; it is possible that the YSO 
happened to form at the correct distance from the exciting star of KR 120 that it now lies on 
the edge of the expanding H II region. On balance though, we prefer the triggered scenario 
as it naturally explains the close spatial association of the H II region and the very young 
Class I YSO. 



5. Conclusions 

1. Our analysis of the HCO+ (J = 3 ->■ 2; 267.6 GHz) and H 13 CO+ (J = 3 ->■ 2; 260.3 
GHz) spectra of the submm/IR source IRAS 01202+6133 shows that the blue-asymmetric 
HCO + line profile is almost certainly due to infall motions. 



2. We have modifie d the analytic radiative transfer models of iMyers et al.l (119961 ) and 



DeVries fc Myersl (120051 ) to allow for fitting of submm/mm line profiles that have both self- 
absorption features and optically thin wings. We applied these models to our HCO+ submm 
spectrum of IRAS 01202+6133 and conclude that it is a young Class I YSO with a substantial 
envelope undergoing slow infall. 

3. Based on its young evolutionary state (t yso < Thii), and its proximity to the H II 
region KR 120 (Sh 2-187), we think that the formation of IRAS 01202+6133 was triggered 
by the expansion of the H II region. However, due to young age of the H II region a " Collect 
and Collapse" scenario is ruled out. 



We would like to thank Dr. Lewis B. G. Knee for his assistance with the processing of 
the JCMT data. 
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1— The H II Region KR 120 (Sh2-187). This Midcourse Space Experiment (MSX; 



Price et al.ll200ll ) 8.3 fim image traces emission from the photodissociation region within the 
molecular cloud surro unding KR 120. Th e red contours correspond to Canadian Galactic 
Plane Survey (CGPS; iTaylor et al.l 120031 1 1420 MHz continuum emission at Tg = 8, 10, 
12 and 14 K. The elongated radio source seen at the bottom of the main H II region is 
extragalactic. The blue contours represent SCUBA 850 /zm continuum with levels of 1.2 x 
10~ 3 , 2 x 10~ 3 , 3 x lO" 3 , 6 x 10~ 3 Jy beam" 1 . Note the location of the Class I YSO IRAS 
01202+6133 on the periphery of the H II region. 
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Fig. 2.— HCO+ spectra of IRAS 01202+6133. The HCO+ (J = 3 ->■ 2; 267.6 GHz) line 
has a classic blue asymmetric line profile. The optically thin H CO + (J = 3 — > 2; 260.3 
GHz) line (with Gaussian fit shown) peaks at the minimum between two peaks of HCO + 
line confirming that blue asymmetric line shape comes from the infall motion. 
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Fig. 3. — Basic model results. The figure shows how the basic Two-Layer model (left) and 
Hill model (right) are unable to fit the extended wings seen in the HCO + spectrum. In these 
plots the systematic velocity of the cloud (Vlsr = —13.7 km s -1 ) has been removed. 



- 14 - 



T, 




Fig. 4. — Modified radiative transfer models. In order to fit the extended wings seen in 
the HCO + spectrum we added a central low optical depth "outflow" region to the exisiting 
two-layer and hill models. The diagrams schematically show the run of temperature with 
optical depth in the 2L-0 model (left) and Hill-0 model (right). 
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Fig. 5. — Best fit 2L-0 model. The model is able to fit the extended wings of the HCO + line, 
but as expected for a two-layer model has some trouble matching the line widths and the 
trough structure. The best fit parameters are given in Table Q] and the systematic velocity 
of the YSO has been removed from this plot. 



- 16 - 



2.0 




-0.5 I , , , , i , , , , i , , , , i , , , , I 

-10 -5 5 10 

V LSR (km s~ 1 ) 

Fig. 6. — Best fit Hill-0 model. The model is able to fit the extended wings of the HCO + 
line and has a better fit to the central portions of the line profile compared to the 2L-0 
model. The best fit parameters are given in Tableland the systematic velocity of the YSO 
has been removed from this plot. 
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Fig. 7. — The molecul ar cloud of KR 120 . The image is 12 CO ( J — 1 — >■ 0) emission from the 
Outer Galaxy Survey (IHeyer et al.lll998l ) integrated between —20.4 < Vlsr < —5.6 km s^ 1 
(contours at 20, 30, 40, 50, and 60 K km s _1 ). The position of KR 120 is shown by the red 
contour at T M2 o = 8 K, and the position of IRAS 01202+6133 is indicated by the cross. 
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Table 1. Best Fit Free Parameters of 2L-0 model 



TO 


ViSii (km s- 1 ) T K (K) ^(kms- 1 ) a NT (km s" 1 ) 


Tout 


(Tout (km s r ) 


11 ± 0.2 


-13.74 ±0.01 9.5 ±0.15 0.02 ± 0.003 0.4 ±0.01 


0.35 ±0.01 


2.3 ±0.03 


Note. - 


Background temperature is set at Tf, = 2.725 K. Uncertainties are 


obtained by 


a Monte Carlo method with 1000 



repetitions. 
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Table 2. Best Fit Free Parameters of Hill-0 model 


Tc 


V LS r (km s- 1 ) ^(kms" 1 ) u(km s" 1 ) T p (K) T (K) tw cw(km s" 1 ) 


6.5 ±0.1 


-13.74 ±0.01 0.07 ±0.012 0.47 ±0.02 7.7 ±0.13 3.2 ±0.22 0.32 ± 0.03 2 ± 0.06 


Note. - 


Background temperature is set at Tf, = 2.725 K. Uncertainties are obtained by a Monte Carlo method with 1000 



repetitions. 



